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The  SAP  IV  Finite  Element  Computer  Program  was  designed  to 
accept  new  elements  into  its  element  library  easily.  The  new 
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TYPE  9  allows  for  material  directions  to  be  arbitrary  for  ease  of 
material  input  descriptions.  The  element  is  modelled  after  the  I 
structure  of  element  TYPE  6;  therefore,  element  TYPE  9  can  1 

degenerate  to  element  TYPE  6. 
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INTRODUCTION 


This  report  describes  research  which  was  directed  at  the 
developnent  of  an  orthotropic  plate  finite  element  for  the 
analysis  of  plate  and  shell-like  structures  which  exhibit 
coupling  between  extension  and  bending.  The  element  is 
especially  useful  in  the  analysis  of  structures  which  are 
fabricated  from  laminated  composite  materials.  The  report  is 
written  so  that  it  would  describe,  for  the  finite  element 
expert,  the  analytical  techniques  utilized  in  the  development 
of  the  element  and  also  would  be  of  use  to  a  program  "user" 
not  having  the  overall  expertise  of  an  expert. 

The  notation  and  the  methods  used  for  the  definition  of 
element  material  properties  have  been  chosen  as  a  result  of 
a  careful  survey  of  the  literature  on  composite  materials. 

The  notation  and  definitions  chosen  are  considered  to  be  in¬ 
dustry  standard  and  are  best  summarized  in  reference  3,  which 
is  rapidly  becoming  a  standard  text  for  the  analysis  of  com¬ 
posite  materials. 

An  industry  standard  computer  program  SAP  IV  was 

selected  as  host  program  to  accept  the  new  composite  plate 
finite  element.  The  SAP  IV  Finite  Element  Computer  Program 
was  designed  to  easily  accept  new  elements  into  its  element 
library.  The  new  element  muot  be  self  contained  since  the 
general  philosophy  and  program  structure  is  "overlayed"  into 
the  computer.  The  laminate  composite  plate  element  is  the  new 
element  to  be  integrated  into  the  element  library.  Called 
TYPE  9,  the  new  elcme.nt  is  similar  to  the  SAP  element  TYPE  6 
in  both  description  and  input.  The  main  difference  is  that 
element  TYPE  9  has  the  ability  to  describe  the  effects  of 
coupling  between  in-plane  extension  and  out-of-plane  bending. 
Element  TYPE  9  is  a  quadrilateral  element  and  is  formulated 
from  quadrilateral  shape  functions  rather  than  from  four 
triangles  as  in  TYPE  6.  Also,  TYPE  9  allows  material  directions 


to  be  arbitrary  for  ease  of  material  input  descriptions.  The 
element  is  modelled  after  the  structure  of  element  TYPE  6; 
therefore  element  TYPE  9  can  degenerate  to  element  TYPE  6. 
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SiJCiION  II 

LJ^-MINATE  COMPOSITE  Fi.AT  PLATE  ELEMENT 


The  laminate  composite  flax  plate  element  is  based  on  thin 
plate  theory  with  the  exclvision  of  transversi  the  ir  deforma¬ 
tions.  Tr.e  following  sections  describe  the  basic  t  irmulat ioii 
of  the  finite  element. 

1.  ELEMENT  POTENTIAL  ENERGY  FUNCTIONAL 

The  principle  of  minimum  potential  energy  furnishes  a 

variational  basis  for  the  direct  rormulation  of  the  element 

stiffness  equations  and  loading  functions.  The  potential 

energy  of  the  element  is  formed  from  tne  sum  of  strain  energy 

(u,  )  and  the  potential  of  all  applied  loads  (V  );  i.e., 

P  D 


The  principle  can  be  stated  as  follows:  Airiong  all  the  displace¬ 
ment  functions  of  adjtissable  form,  those  that  satisfy  the 
element  equilibrium  conditions  ir.ake  the  potential  energy  func¬ 
tional  obtain  a  stationary  value.  Thus, 


£  .1  =  6u  +  6V  -  o  ( 1 ) 

P  P  P 

where  6  is  the  : irst  variational  operator. 

It  can  be  shown  that 


id 


p 


6e  dV 


(e) 


where  o'"  is  a  vector  of  stress  components, 

e  the  corresponding  vector  of  strain  components  and 
V  the  volume  of  the  element. 

Note:  All  vectors  will  be  underscored  with  a  straight  bar  and 

matrices  will  be  underscored  with  a  tilda.  The  super- 
T 

script  of  the  vectors  and  mati-'ices  designates  the 
matrix  is  transposed. 


The  cox'responding  first  variation  of  the  potential  forces 
becomes 

5Vp  =  -  /"  6u  dV  -  f  ^  6u  dS  (4) 

V  S 


where  b 

i-^epre 

sents  the 

element 

body  forces , 

IS  a 

vector  or 

surface 

tractions  applied  on 

Sux^idcc  S,  dnei 

u 

is  a 

vector  of 

element 

displacements . 

here  that  the  surface  traction  integral  can  be  used  to  include 
the  point  concentrated  forces  on  the  boundaries  of  the  element. 

The  elements  of  the  strain  potentials  of  equation  (3)  will 
eventually  lead  to  the  element  stiffness  and  initial  load  vec- 
■^ors  and  the  elements  of  the  applied  load  potential  of  equation 
(^)  will  produce  the  various  element  vectors. 

2.  QUADRILATERAL  SHAPE  FUNCTIONS 

The  element  formulation  is  a  geometrically  linear  quadri¬ 
lateral  containing  the  four  corner  nodes  as  shown  in  FIGURE  1, 


FIGURE  1  -  Quadrilateral  Element  Geometry 
a.  Geometric  Shape  Functions 

The  element  shown  in  FIGURE  1  is  described  in  the  local 
coordinate  of  the  element  and  all  material  reference  is  made 
with  respect  to  the  element  local  x  axis.  The  element  area 
domain  can  be  described  by  using  a  polynomial  as 


4 


(5) 


T 

X  =  ^  B 
T 

y  =  i  B 

where  x  and  y  are  the  local  coordinates  as  in  the  element 
domain , 


T 

d"  -  Ll,  r,  s,  rs j  (6  > 


the  row  vector  of  polynomial  coefficients  a  vector  of 
generalized  coefficients,  and  r,s,  the  element  natural 
coordinates . 

The  generalized  coefficients  can  be  solved  for  by  evaluating 
the  polynomials  at  the  vertices  of  the  element.  Therefore, 


X 

y 


X 


T 

H  Z 


(7) 


T 

where  H  contains  the  terras  of  the  shape  function  of  the 
element;  x  and  y  are  vectors  containing  the  element  ver- 


x'  =  lx^,  x^,  Xg,  x^^J  (8) 

T 

y  =  Ly^,  y2,  y3,  y4J 

The  terms  of  the  shape  function  can  be  described  by 

h^  =  ^  (1  +  r^r)  (1  +  s^s)  (9) 


where 


r^'^  =  L-1,  1,  1,  -IJ  (10) 

=  L-1,  -1,  1,  IJ 

define  the  natural  coordinates  ot  the  elements. 
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The 

can 


mappinji  of  the  element  geometry  and  displacement  functions 
Le  obtained  by  defining  the  Jacobian  transformation  as 


3  =  J  3 

-n  -  -I 


(11) 


where  J  is  the  Jacobian  matrix  defined  as 


and 


a 

-n 


3X 

a 

ay 


(12) 


(13) 


are  the  first  derivative  operators  in  the  natural  and  local 
reference  frames  respectively. 

i.'ote:  The  " , subscript  implies  "partial  differentiation  with 
respect  to".  The  inverse  transformation  is  obtained  by 

=  j.  5  '1 


where  J*  is  the  determinant  of  the  Jacobian  matrix  given  as 


J* 


(16) 


and , 


6 


It  will  be  necessary  to  obtain  second  derivates  in  the 
local  reference;  therefore 


E  In 

—  — n 


+  F  ^ 


(1 


represents  the  second  partial  operator  in  the  local  reference 
given  as 


3xay 


and  the  natural  set  as 


3 


2 


-n 


/ 


3r3s 


(1 


(2 


The  E  matrix  is  defined  as 


where 


e  .  . 

-i: 


1  -  T  d*  I 

J*  -n  -i 


(22) 


with  gt  being  the  ith  row  partition  out  of  the  G  matrix  and 


3 

-n 


9n  -  7*  ^n  ^ 


(23) 


The  F  matrix  is  defined  as 


hi 

f  ^ 
-22 


f  ^ 

-12  -• 


(24) 


where 


with 


^  T 

-ij 


—  2 

m 

ii  ij 

(25) 

hi 

h2  ° 

(26) 

o 

®il  ®i2 

- 

>f  the 

G  matrix. 

b.  In-Plane  Displacement  Shape  Function 

The  plate  element  is  assumed  to  have  in-plane  deformations; 
therefore  the  variation  of  x  and  y  displacements,  u  and  v  res¬ 
pectively,  can  be  expressed  using  the  same  shape  functions  for 
the  geometry  as  in  the  previous  section.  Then 


h"  % 
s’'  % 


(27) 
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are  the  domain  displacements  of  the  element 


where 

T 

3u  - 

m 

•■“l  “2  “3 

(28) 

and 

% 

“2  "a 

contain  the  in-plane  model  displacements.  Therefore,  it  is 
assumed  that  in-plane  displacements  vary  linearly  within  the 
element . 

c.  Transverse  Displacement  Shape  Function 
The  plate  element  defined  by  thin  plate  theory  must  have 
a  transverse  shape  function  to  allow  for  proper  bending. 
Therefore  it  is  assumed  that  the  shape  function  polynomial  is 

T  I  2  232  233  3 

<J>  =  Ll,r,s,r  ,rs,s  ,r  ,r  s,rs  ,s  ,r  s,rs  J 

(29) 

The  natural  degrees  of  freedom  allowed  per  node  for  bending  are 

aoi  =  '"’r  ”’s-'  i 


Letting 


3o  “ 

where  is  the  full  set  of  degrees  of  freedom  in  the 

natural  reference  of  the  element,  ip  is  the  matrix 
$  evaluated  at  the  nodes  and  defined  in  reference 
2  and  6  is  the  set  of  generalized  nodal  coefficients. 
Then  the  transverse  displacement  w  becomes 

”  =  ^n  ^0  " 

where  ^  is  the  inverse  of  and 

II  i  : .  I  I  ii  I  r'.irr  ;ve  r’: ;  h.  i  pe  fu  tn  ■  I  i  <  ih  ,  both  o  I  wti  i  ch 
,  I  r  ■!  .  S  I  I  1 1 ' I  .11  t  •(  •  I  <  ■  t  ■  1 1<  ■< '  2  . 
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The  transverse  displacement  shape  functions  are  defined  in  the 
natural  reference  of  the  element  to  allow  for  ease  of  develop¬ 
ment  since  second  derivatives  must  be  taken.  Since  a  global 
transformation  by  node  is  to  be  performed  at  a  later  stage, 
the  transformation  by  node  from  the  natural  to  local  coordi¬ 
nate  can  be  made. 

3.  PLATE  STRAIN  FUNCTIONS 

The  classical  assumptions  of  linear  thin  plate  theory  are 
made,  essentially  reducing  the  three-dimensional  equations  of 
tjlacticity  to  a  two-dimensional  set  of  plane  stress  equations, 
for  the  elastic  continuum  of  the  plate,  the  following  assump¬ 
tions  are  made: 

•The  thickness  (h)  is  small  compared  to  the  dimensions  of 
the  plate  in  the  x  and  y  directions. 

•A  line  element  through  the  thickness  remains  normal  to  the 
mid-plane  surface  under  all  states  of  deformation,  inde¬ 
pendent  of  its  translation  or  rotation. 

•The  plate  can  be  isotropic,  orthotropic  or  comprised  of  a 
number  of  orthotropic  laminae,  where  each  lamina  obeys 
Hooke's  law. 

•The  displacements  u,  v,  and  w  in  the  x,  y,  and  z  direc¬ 
tions  respectively,  are  small  when  compared  to  the  plate 
thickness . 

•The  reference  axis  is  taken  as  the  middle  of  the  plate  at 
h/2,  h  being  the  total  plate  thickness. 

•The  normal  strain  in  the  z-direction  is  assumed  to  be 
zero,  giving 


e 


z 


w . 


0; 


therefore,  the  lateral  deflection  is  given  by. 


w  =  w  ( X ,  y )  . 
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•St.  Venant's  principle  applies.  That  is,  local  deformation 
occurs  in  the  area  of  applied  loads  while  at  distances 
away  from  the  load,  the  deformation  state  is  not  grossly 
affected . 

•Transverse  shear  deformations  are  neglected; 

Y  =  Y  =  0  . 
xz  yz 

•Displacements  are  linear  such  that 

u  ^  UQ(x,y)  -  zw,^  and  v  =  VQi.x,y)  -  zw,^, 

whe.-'e  w,  -  -6  w,  -o  and  and  are  the  in-plane 
’x  y’  y  X  0  0 

displacements  of  the  middle  surface.  The  rotations  about 

the  X  and  y  axes  are  given  by  6  and  Q  ,  respectively. 

X  y 

a.  Midplane  Strain  and  Curvatures 

The  mechanical  strains  associated  with  plate  stretching 
nd  bending  can  be  written  as 

e=e+z<  (5j) 

-  -o  - 

where  the  mid-plane  strains  are 

(  34  ) 


the  plate  curvatures  are 

XX  j 

(35) 

yy  ( 

’  xy 


1  L 


with  u  and  v  being  the  in-plane  displacements  and  w,  the  trans 
verse  displacements.  The  thermal  strains  can  be  written  as 


c.  =  ST  +  2aT 
-t  -  o  -  g 


(36 


where  d  is  a  vector  of  thermal  expansion  coefficients 
relative  to  mid-plane  strains,  the  element  mean  temp- 


erence  and  the  element  thermal  gradient 


eraturc  an 
through  the  plate  thickness. 


b.  Strain  Displacement  Functions 

The  connection  between  strain  and  displacement  is  made 
realising  that  the  in-plane  and  transverse  displacements  have 
been  made  relative  to  a  set  of  nodal  displacements.  Equation 
(33)  can  be  written  as 


^  2i  +  zBq 


(37 


where  is  the  in-plane  strain-displacements  relative  to 
(the  in-plane  nodal  displacements); 

Bp  is  the  transverse  strain-displacements  relative  to 
(the  natural  transverse  nodal  displacements). 

The  in-plane  aisplacements 


3i 


(38 


and  the  becomes 


wher  e 


(40) 


The  out  of  plane  displacements  are 


i  3o4  j 

I  j 

where  the  sub-elements  of  the  partition  are  defined  by 
equation  (30)  and  Bq  becomes 


(41) 


h  ?n  *  r  iJn’ 


(^2) 


where 


r-: 


0 

0  -1 

0  0 


0 

0 

-2 


(43) 


-n 


3  R- 

-n  - 


y-T 


(-H  ) 


ana 


-T 

H 

-nn 


2  — T 

2n  S 


(4b) 


Mote:  In  equation  (37)  the  z  variable, which  is  the  plate's 

normal  coordinate,  is  maintained  distinctly  since  it  is 
independent  of  the  in-plane  variables.  Later,  when  the 
strain  energy  is  formed,  the  z  variable  will  integrate 
through  the  thickness  and  merge  into  material  property- 
matrices  . 
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;''L.-\NE  STRESS  COMPONENTS 


The  stress  components  for  a  thin  plate  can  be  written  in 
vector  form  as 

0  =  C  (e  -  e^)  (46) 

where  C  is  'che  material  matrix  described  in  the  plate  local 
axes  and  is  expressed  as 

C  C  R  (47) 

'c  '■  ~e 

witii  C  being  the  material  matrix  in  the  principal  material 
directions  of  the  fibers  and  being  the  strain  trans¬ 
formation  matrix  from  element  focal  coordinates  to  prin¬ 
cipal  fiber  directions.  The  elements  of  equation  (47) 
are  found  in  reference  2. 

The  elements  of  the  thermal  strain  involving  the  thermal  coef¬ 
ficients  are  defined  as 


i  a  . 

“e  ~ 


(48) 


Once  the  material  matrix  is  defined,  the  elements  of  the  material 
...lasticity  matrix  can  be  defined  as  in  the  following  section. 

1.  MATERIAL  ELASTICITY  MATRICES 

Thv  material  coefficients  are  defined  with  the  use  of 
equdtion  (3)  in  a  slightly  different  form: 


e  dV 


(h9) 


Since  the  local  z  dimension  is  small  compared  to  the  x  and  y 
plate  dimensions,  it  is  convenient  to  define  the  stress  resul¬ 
tants  and  moment  resultants  as 


A 


adz 


(SO) 


1  4 


ina 


/  2 

t 


z  d  z 


rhen,  a  new  stress-strain  matrix  can  be  defined  as 

-  i-  _  1-0 

'  L  ut  pj(. 


where 


=  /J 
=  /  ? 


C  dz 


C  zdz 


C  z  dz 


Letting 


f  A  B 


D 


-T  5t 


equation  (49)  can  be  written  as 


2’'  £  dv  =  /^  -m  ^  <1A  -  /a  5 


-  -Tl,Tg/<iA 
(59) 


L 


where 


At 

-  /, 

1  O  ) 

a 

dz 

(60) 

-T 

-  L 

C 

a 

zdz 

2 

(61) 

PI 

A 

c 

a 

z  dz  . 

(62) 

The  material  A,  B,  D  matrices  and  the  thermal  load  coef¬ 
ficients  A^,  B  and  D  can  be  related  to  laminar  material  by  posi- 
tion  t  in  the  material  build-up  as 


A  = 

L 

z 

i  =  l 

c.  (t.-t.  ,) 

~i  1  1-1 

(63) 

B  = 

1/2 

Z,  A  - 

(64) 

D  = 

1/3 

^  -  1 

(65) 

A„  = 

-  i 

L 

z 

i  =  l 

Si  “i  'A  -  A 

(66) 

-T 

1/2 

y  C-  a.  (t^ 

-  'li> 

(67) 

-T 

1/3 

y  C.  a.  (t? 
i=l  ^ 

-  ^Li> 

(68) 

where  the  subscript  "i"  implies  coefficient  evaluation  at 
laminae  level  "i"  and  L  is  the  total  number  of  fiber  lamina 
levels . 

6.  COORDINATE  TRANSFORMATIONS 

The  element  information  is  initially  determined  in  the 
natural  coordinates  of  the  plate  since  it  is  quite  easy  to  expres 
all  loading  and  stiffness  information  in  that  reference.  Ulti¬ 
mately  the  information  must  be  transformed  to  local  coordinates 
(x,y,z)  and  also  to  global  coordinates  (X,Y,Z).  The  following 
sections  describe  the  transformations. 


-i .  Natural  tr-  i^O‘-'ai  Tra/iu:  ormat  ion 

The  iidturai  coordinate  variables  per  node  are  defined  as 

/  \ 


J  "  '■ 

'■^ni  v  /  I 

I  I 

wnere  is  the  rotational  degree  of  freedom  normal  to 
plate  at  node  "i". 

he  transformation  matrix  required  becomes 


^ndi 


where 


t-  ,  ■  =  0 

-nil 


I 

L  0^  0^  1_ 


LuVWU  0  GJ 
X  y  z 


with 


i:i- I 


.Voting  that  the  rotation  degrees  of  freedom  are  definec  as 


-^’x  I 


ihv;  r  raatrix  become  a 


Q-: 


(74) 


The  complete  natural  to  local  transformation  therefore  becomes 


~  n£l 


~n£2 


-n£  J 


■n£4 


(75) 


b.  Local  to  Global  Transformations 

The  local  to  global  transformation  quantities  are  somewhat 
more  difficult  to  obtain  since  the  transformation  involves  the 
local  coordinates  of  a  quadrilateral  element.  Obviously  only  3 
points  define  a  plane;  the  fourth  point  of  the  quadrilateral  is 
unnecessary.  The  fourth  point,  however,  may  not  lie  in  the  same 
plane  as  the  other  three  points.  Therefore  local  transformations 
by  node  are  determined  and  are  averaged  to  obtain  a  general 
transformation  used  in  determining  the  element  coordinates  and 
in  transforming  element  matrices  where  applicable. 

Defining  the  nodes  of  Figure  1  as  i,  j,  k  and  1  and  allowing 
this  sequence  to  permute,  the  element  normal  coordinate  at  node 
"i”  which  also  permutes,  is 


zi 


V .  -  X  y . , 

I21 _ Zll 

IV  .  .  X  V.  J 
-:i  -ii 


(76) 


where  the  "X"  symbols  denote  a  cross  product  of  two  vectors, 
and  V . .  implies 

-31 
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z.  -  z.\ 


V, itr.  X,  Y  and  Z  being  global  coordinates. 

A  global  material  reference  is  defined  as  a  global  vector 
defined  as  input  which  locates  the  general  local,  "x"  of  all 
elements  referred  to  that  vector.  All  material  properties  are 
defined  relative  to  this  x  coordinate  which  becomes  the  ele¬ 
ment's  local  X  coordinate.  The  element's  local  y  coordinate 
at  each  node  can  be  calculated  as 


e  .  X  V 
zi  -m 

6  -  X  V“I 
zt  ~m 


(7c) 


’inally,  the  local  x  coordinate  at  each  node  is  determined  as 

(79) 


g  .  X  ^  . 


XI 


le  .  X  §  . 
f  yi  zii 


The  local  to  global  transformation  at  node  "i"  becomes 


X 


i 


(80 


Le; 


where  the  direction  cosines  of  each  coordinate  are  placed 
in  row  order  in  the  transformation  matrix. 

The  average  transformation  used  to  determine  the  local 
coordinates  is  obtained  by  first  averaging  the  nodal  normals  as 


e 


z 


(81) 


and  substituting  into  equations  (78,79  and  80)  to  produce  an 
average  t^^^. 
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"he  global  degrees  of  freedom  can  be  defined  by  node  as 

V 

■gi 

e, 

e. 


where  U,  V  and  W  correspond  to  global  displacements  rela¬ 
tive  to  X,  Y  and  Z  respectively  and  6*s  corresponds  to 
global  rotations  about  X,  Y  and  Z  respectively. 

The  complete  transformation  becomes 


where 


~lgl 


■ig2 


*g4 


(83) 


.igx 


(84) 


and  q„  and  q  are  the  complete  list  of  degrees  of  freedom 

JC  -“g 

per  element. 

If  the  quadrilateral  element  is  perfectly  flat  in  its  space, 
then  the  t„  •  becomes  exactly  t„  .  If  it  is  not,  then  the 

,.£gi 

element  space  appears  as  a  curved  space.  This  effect  should 
allow  the  element  to  behave  as  a  shallow  shell. 


c.  Natural  to  Global  Transformation 

The  element  stiffness  will  be  transformed  from  natural  to 
global  coordinates  directly.  Therefore  that  transformation 
becomes 
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2 


n 


2g 


(85) 


where 


~ng 


(86) 


Since  both  original  transformation  matrices  are  partitioned 
diagonally,  the  nodal  transformation  matrix  is 


t  .  =  t  „  .  t .  . 

~ngi  ~n£i  -.£gi 


(87) 


which  contains  many  off-diagonal  zeroes.  Therefore  it  is  con¬ 
venient  to  define  transformations  by  displacement  and  rotation 
degrees  of  freedom.  That  is,  the  displacement  transformation 
at  node  i  is 


^Di 


!i 


(68) 


and  the  rotation  transformation  at  node  i  is 


This  modification  will  save  transformation  operations  later  on. 

7.  ELEMENT  STIFFNESS  MATRIX 

The  element  sciffness  is  easily  defined  in  the  natural  coor¬ 
dinates  of  the  plate,  given  the  strain  displacement  functions. 
Once  this  stiffness  is  determined,  it  can  be  augmented  with  a 
scaffolding  or  artificial  torsional  stiffness  for  the  plate's 
normal  degrees  of  freedom.  Finally,  this  stiffness  can  be 
transformed  to  local  and  then  to  global  coordinates  for  assembly 
into  a  master  stiffness  matrix.  The  following  details  the  above. 


21 


a.  Plate  Element  Stiffness  in  Natural  Coordinates 
Defining  a  new  set  of  degrees  of  freedom  as 

the  strain  components,  as  defined  by  equation  (58)  become 


B 


~n 


§ 


n 


(91) 


where 


(92) 


The  second  integral  of  equation  (59)  can  be  used  to  define  the 
element  stiffness  in  local  coordinates  as 


E„  e  dA 
~m  - 


K  q 
~n 


where 


(93) 


R 


~n 


E 


~m 


B  dA 
-n 


(94) 


For  convenience  of  computation,  the  material  matrix  E^  in  equa¬ 
tion  (94)  is  Cholesky  factored  as 


T 

E  =  U  U 
~m  -  ^ 


(95) 


where  U  is  a  upper  tri-diagonal  factoring  matrix. 

This  allows  equation  (94)  to  be  written  in  a  more  efficient 
form  as 

K  -  f  (U  B  )’’’  (U  B  )  dA  (96) 

~n  J  ~  ~n  ~  ~n 

A 

which  allows  the  triple  matrix  product  to  be  replaced  by  a  sim¬ 
pler  transpose  symmetric  product.  This  process  is  especially 


2? 


efficient  since  the  Cholesky  factoring  is  performed  (at  most) 
once  per  element.  Also,  numerical  integration  is  to  be  per¬ 
formed.  Savings  will  occur  at  each  integration  point  after 
the  first. 

b.  Artificidl  Torsional  Stiffness 

The  flat  plate  theory  does  not  have  any  mechanism  to  di¬ 
rectly  include  twisting  of  the  plate  normal  to  the  plate  sur¬ 
face.  Therefore,  if  two  coplanar  elements  are  assembled  at  a 
comimon  node,  a  .-singular  stiffness  exists.  To  avoid  this,  an 
artificial  or  scaffolding  stiffness  is  added  to  the  normal 
rotational  degree  of  freedom  6^.  There  is  no  change  in  the 
system  equi libriu.m .  For  convenience,  this  is  performed  at  all 
nodes  of  the  element  since  it  would  be  difficult  to  determine 
coplanar  effects  in  general.  This  does  change  the  overall 
element  equixibrium.  If  the  amount  of  artificial  stiffness  is 
Kept  small  and  the  local  rotational  stiffness  effects  are  in 
equilibrium,  then  the  error  can  be  minimized.  Defining  a  vec¬ 
tor  of  normal  rotations  at  the  nodes  as 


the  artificial  torsional 


■fness  matrix  relativt 


to  d  bu.' comes 


where  f  is  an  input  scaling  factor  wr.ich  can  vary  an 


0  <  f  <  1 


(99) 


and  can  be  set  in  l.E-8  increments,  C  is  an  artificial 
coefficient  estimated  from  element  bending  stiffness 
coefficients  and  element  area;  i.e., 

C  =  MIN  (D(l,l),  D(2,2))  *  AREA  (100) 

with  the  D's  defined  in  equation  (55). 
c.  Natural  Stiffness  Matrix 

The  degrees  of  freedom  q  and  6^  defined  by  equations  (90) 
and  (97)  can  be  merged  to  degrees  of  freedom  described  in 
equation  (69).  This  requires  the  re-ordering  of  stiffness  coef¬ 
ficients  of  equations  (96)  and  (98)  to  produce  a  natural  stiff¬ 
ness  matrix: 


K 


~n 


MERGING 

REORDERING 


(101) 


relative  to  q^. 
d.  Global  Stiffness  Matrix 

The  global  stiffness  matrix  K  is  formed  by  transforming 
from  natural  coordinates  to  global  using  equation  (85).  The 
transformation  is  formed  using  equation  (1)  realizing  that  the 
strain  energy  is  invariant  relative  to  any  coordinate 
reference.  Therefore, 


U 


P 


Using  equation  (85)  produces 


K 

~r, 


(102) 


(10  3) 
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inefficient,  gs[h:(' i  ,i  1 1  y  since  fjj,  hii'.hly  (jiay.onul.  !,f  :  .cicncy 
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Can  n&wt;\'>^r  dc  eiitrcita  -jV  part  it  ionintr  k  and  K  into  J  /. 

~  n  ~  p 

n  P-  ^ 

s'cP-.natrices  laoellea  ^  and  K^-v"  where  i  and  i  range  from 


to 


rhen 


-  J 


T  -  K.  .  r  . 


(104  ) 


wiiere  t^-  matrix  relates  to  liquation  (88)  when  i  eqaal 

1,  3,  5,  7  ana  relates  to  equation  (89)  when  i  equals 

2,  4,  6,  3. 

Additional  etiicieriCy  is  oOtaincc  wnen  tnu  tripo-c  rniati'ix  prcauci 
IS  performed  such  that  the  right  portion  matrix  multiply  is 
first  formed  ana  positioned  back  into  K.*?.  Then,  the  left 
multiply  is  formed  with  the  resulting  product  ana  placed  into 


?.  Numerical  ixteokaticx  of  area  functions 

The  elements  in  the  matrix  of  equations  (Sh)  and  (96)  are 
very  difficult  to  integrate  exactly,  therefore  approximate 
numerical  integration  can  be  performed  with  sufficient  accuracy 
for  i ustif ication .  Gauss-Legendre  Numerical  Quadrature  has 
seen  selected  to  perform  the  integration  o:  the  stiffness  coef¬ 
ficients  as  well  as  other  area  functions.  Equation  (b4)  can  he 
rt'u t' i 1 1 e n  and  transformed  relative  to  variables;  and  limdts  of 

.  *  L  o  r  a  L  -I  vv  r.  '..j  c:  s 


y  y  (r,s)  E  (r,s)  J-  cm 

-i  ~n  ~m~n 


(  .  .  i.  ) 


An  element  of  thi.,;  matnx  can  be  written  as 

fl  rl 

k.  .  =  J J,  ;  .  .  (r,s)  ;n"ds 
1  '  -  i  - 1  1  ] 


1. 1  6  c  ) 


where 


f:  •(r,s)  -  E  (k,i)  I  i:  (K,n  :■  (^,i) 


-n 


(.87) 


I 


The  stiffness  coefficient  can  then  be  approximated  as 

hi  \  <1 

where  n, ,  n„  are  the  number  of  Legrendre  root  evaluation 

X  i 

points  in  the  r,s  directions  respectively, 
r.,  and  s^  ai'e  the  roots  of  the  Legrendre  polynomial, 
and  are  the  appropriate  Gauss  weighting  factors. 


:  or 


r  or 


^  ’VT 

Wk  =1,1 


(109) 


n  =  3  , 


Wk 


Ys”  Vi~ 

”  ,  0,  5 

5  8  5 

9  ,  9,  9 


(110) 


9.  ELEMENT  MASS  MATRIX 

A  consistent  mass  matrix  relative  to  the  in-plane  variables 
in  one  coordinate  can  be  written  as 


p  t  H  H  dA 


(111) 


and  a  corresponding  lumped  mass  matrix  can  be  formed  by  summing 
the  rows  of  the  consistent  mass  matrix  as 


M 


I 


p  t  (  dA 

‘^A  ■ 


assuming  t  and  p  constant  and  realizing 


(112) 


26 


Tht:  c-j",v.Gner.t L  g  -ir’e  applied  to  the  trans latoi-'y  dugree^  c.f 
.  rtedc^.  Ill  ill  dire  cions,  per  node.  A  rotary  inertia 

eiiect  can  be  includea  by  an  approximation; 


:n<,-  piGDCix  mumped  mass  vector  can  ce  formec: 


(113) 


1 


I  :  ‘ 


(114) 


wl'icre  i  IS  a  G  1  unit  vector. 

1C.  ELEMENT  LOAD  VECTORS 

The  element  load  vectors  are  established  from  element  prop¬ 
erties  such  as  material  constants,  temperatures,  pressure,  n.aGs, 
area  and  acceleration  constants.  The  following  sections  des¬ 
cribe  the  load  vectors  ‘developed . 

11.  Tiier'ma)  1  gic  Vector 

Tne  there  intcp.ral  of  equation  (59)  c:an  ta-  usee  t^  def.r.'.- 
t'ne  element  t .-.(..r.m.il  load  vector  using  equcition  (Si): 


“  =  a; 


dA 


^  ro( 


There  f  ore , 


(110) 


is  tfic  therm. il  vei  tor  relative  to  tiie  n..it'i!.i]  coordinatc-s 


a,,,  I:  I  S' 


ri:::,  v. 


iiM  '■  I  in  .  I  t  I  I  ij  .  : 


7 


(117) 


Ir  =  Et  {t  j  • 

The  global  thermal  vector  can  be  determineu  by  applying  the 
natural  to  global  tranuf ormation , 


T  "  F" 
~ng  -n 


(  i  1 0  ) 


b.  PreL'. sure  aoad  Vector 

The  pressure  load,  by  node,  is  formed  for  thi^  shape  terms 
associated  only  with  the  geometry  and  is  appliec  in  the  z  (plate 
normal)  direction.  Therefore, 


f  p  K  dA 
u'A 


(119) 


A  pressure  loaa  vector  can  be  formed  relative  to  local  coordi¬ 
nates  as 


(120) 


The  global  pressure  load  can  be  formed  as 


^  ^  pP' 

.£g  ££ 


(121 ) 


which  transforms  the  normal  traction  into  a  gj-obal  traction, 
c.  Constant  Acceleration  Load  Vector 

The  acceleration  vector  can  be  computed  from  the  ma^■s  vu  cto: 
Je fined  in  equation  (114);  i.e.. 


(  1  2  2  ) 


SECTION  :ii 

MODxC.CA'CTONS  TO  SAP  IV  COMPUTER  FROGRAI' 

Tht.'  SAr  TV  cc  .'nruter  program  w.^a  modified  to  accept  the 
i^'ornpod :  t<-  eii.mer.t.  The  following  cectiona  describe  changes  to 
i;a--  exxsling  pe-'ogr'am  as  well  as  new  routines',  of  the  composite 
■p.  Ic.  ft  e  Lenient. 

SAi'  TV  STRUCTUdU 

The  main  onanpes  to  the  SAP  IV  program  occur  in  tne  elemen 
_icrary  centre!.  They  occur  in  the  olem.ent  generation  portion 
of  the  program.  Figures  2  and  3  depict  the  major  routines  used 
.y  SAP  to  perform,  a  static  analysis  of  its  elements.  The  EFT'/;’ 
j  '-..tine  call'd  a  new  routine  called  Pa.ATL,  wfiichi  is  usee  to  ca.  l 
(VcKL/\Y  (8).  (stress  recovery)  Once  OVERLAY  (6)  is  called,  a 
call  i.s  made  to  CFLATE,  the  main  routine  of  the  composite  -plate 
fi.iite  element.  Figures  4,  5,  6,  and  7  describe  the  flow  of 
I'outines  used  by  CPLATE. 


^  I 


I 


DE  SCRIPT  IT' A 


Input  title  ana  master  control  cdra 
Input  nodes 

Write  8  ID,  X,  V,  Z,  T 

Cali  e lenient  type  (new  composite 
plate  element  called  here) 

WRITE  1  Element  stress  matrices 
for  recovery 

WRITE  2  Element  stiffness  matrices, 
mass  and  load  vectora. 

Input  Concentrated  ' oads/masscs 
WRITE  3  Concentrated  load;-;/ 
masses  for  assembly 

READ  2  Read  stiffness  roi'  asocmblv 
WRITE  4  Assembled  ^^tiffr.ess  and 
l.oads  in  blocks 
WRITE  3  Assembled  Mass  Vector 

Call  appropriate  solution  type 
0  Static 

1  Nodal  Extraction 

2  Force  Response  Analysis 

3  Response  Spectrum  Analysis 

4  Direct  Integration 


BASIC  PROGRAM  EI.OW 


I;  Kx  !' 1  H  i  . 

Static  ‘_,oluti:.;i  option 

Solution  or  o iopiacc-ment  equation 
READ  4  Read  -'tiff ness  and  ioac„ 
decomposition 

WRIT;!  2  Solve  for  displctceme.-.t s 

READ  8  Read  ID  identif ication 
n.atrix 

READ  2  Recover'  displacements 
Print  displacements 

READ  1  Element  stress  tarricGS 
Stress  recovery  and  Printing 


C  SOLUTION 


AND  RECOVERY  .STAGE 


DESCRIP'';:ON 


Composite  plate  routine 

Input  material  Properties  in 
tables 

Input  f.lobjl  material  vector.s 

Input  loan  factors  TLO  as  a 
■unction  of  thermal,  pressure 
and  constant  acceleration 

Input  element  definition 

Form  ABD  matrices  and  calcu¬ 
late  0^  thermal  stress  vector. 

Kodify'^E  =  Uf  U 
~m  ~  ~ 

-Petdrmine  local  to  global  trans¬ 
formation  and  locax  element 
coordinates . 

FROSN  -  forms  vector  cross  product 
ST.RAhS  -  special  transpose 
KGI.  -  used  to  transform  global  to 
local  coordinates 

i 

Set  element  K  ,  B  ,  A,  M  matrices 
_ _  -n’  .n’  -g 


FlClld;  U  FLOW  DIAGRAM  FOR 


Ui'tJ  ^ 


~  .L  y  l\'- 

=  i  ,  NO 


r  i.oiv 


I 

j.oop  on  :'-o  T'oints)  in 

^ '  a  C  i'l  <11 1  PC  C  "lC:  o  n 

Ijerermlne  natcr-al  .nape  funct  -cn 

(II,  I?) 

-n 

Accumulate  B  over  all  Gauss 
.  ^  ~  n 

points 

'IKuLT  -  forms  ..pecial  multiply 
of  L'  B.- 

ATAx  -  form  special  nyntnetric 
multi:  Iv  (G  B„)'(L'  fj-j) 


i  '■  1 1  ’  iTi  o  v.  p  i  1  a  e 

mat 

rix 

prCvICCt  taHC 

..tores  result 

i  n 

K 

r. 

Lower  syrrjnotr' 
form^  K  (11, 

12) 

n  <  ■  V ; 

4.  V.  ’  r '  m  u  j.  t  i.  p  -i  ’  / 

w(ii)'-w(i:')*K 

Lower  syrrimetr 

ic 

a  a  c 

-  ads  intr.- 

.•ration  resul 

to 

accumulat'S  r 

foi-  local  Clement  stiffness. 

.Betermr.es  area  ciftribution  vec¬ 
tor  A  by  nocn. 


lAGRA'.M  FOR  Gi''l,A'.'F  (G'  u'f’) 


L21.IJ 


;pi:-c:..!i  mair.x  rnui’eipi',' 
-m 

i’oT’m  trlen'toiii  .  mat 

relat  Lve  ta  /, i>  ■  .  'lii  r' 

a  L  i.ocar*‘ari 

.’M.mbers  foi’  element  fio 
connectivity. 

a  j.  c  u  1  a  t  e  banc  Vv  i  o  t  'n  c  >  i 
n  T  i.  1 1  n  c  a  b  ,  x  o  a 
anc  masc  vecto 
as  aem.bly 

•jtcess  matrix 
recovery 

c on T  r o  1  t  o  ma  i 


iIA>lRAM  }'0k  CI  ;.AV:  (C''!:'?) 


WRITi'  1' 


WRITE  1 


I 


..  -  .r.'.ir.,.  .A  deA:cr’ i  pt  ion  o.  the  routines  usc-a 

. ..  ti  '.'Late  overiay. 

;:.a  ^ea  by  CFLATI.'  to  fill  tfiO  A,  n  and  ’0  o: 
-:o,  ■  and  55  material  .tiatriceu  r  i-om  the  input 
■  lunL.A  arrayn.  Cqaationu  CO,  C'i  ,  and  G2  .^re 

.  -  :  Tin  thernu;’.  lo.id  vectors,  and  thermal 

equation  r  a  7 . 

Cdlleii  by  L'K.A;’  to  j’acto-  cine  matervic:! 
into  an  upper  triunsular  matrix  as  describea  i:i 


-  .  iia'  sut  ill-,  called  by  Cr'LfiTii  to  deter.mine  the  local  to 

.  ;  S' ■  1  :t";nation  matrix  as  described  oy  ecu.itions  76 

.  .  i'hfc-  routine  cheexs  for  proper  area  definitions, 

.  -  :  ..  ''Am  i::',  called  oy  i,0CAL  to  perform  vector  cj-tsis 

:  .  in  ■.-.nation:;  70,  78  and  79.  I'he  x^esulting 

.  .'t  r  poMe  r.  ,  are  norma  li;',od  to  unit  vectors. 

-  ..  ."u.u  ou  t  iia'  Called  by  i.OCAi.,  to  perform  an  in-place 
:n...;e  :.,.o.rix  r-aAspose  •.sf  the  local  to  global  transfor- 
",  ao:.  •  -N.  i  on  oO  . 

-  us;  ',;'ir;e  ..led  by  L,0d.‘i,  so  trannfcnrm  glob-il  element 

:  A-:  i-  a'i  element  -jAUj's-i  inates .  Tra.-isformai  i  on 

r  j";  ix  e  . I  t.-, by  equ-itioM  80. 

-  ...  _  T.  i.:,(  -  All '.  v  to  aCt  matrix  arrav  space  to 


n  \  V  .  J  .  ( .  a  :  . 


'  :  icai  iy,  it  ir,  ..sed  to  set  .matrix  :sp.;iCt' 


-  .a  :  outine  caile-d  by  Crq.Aidi  to  ;d;rm  composii'’"  plaie 
-  1-.;, train  u  i.tnla  -eiiient  tunction  :■  aS  dcscrilttd  hv 

-  :  i 

■  u,;t,  :o:.  ni  .in  l  u..es  ecju.ition:.  ''  ii'.'’  .u,'h  75. 

;  -  il;  :’iL.led  bv  Ci'UAl:’  t,.  ,nAjm;jlat>'  1:a'  r-  :  a  i  i  l 

- ' ri  at'iA  r  .  ...  t:i  ,  ; .  ;(1  r -A  t :  j  r  t-  -..A  al  ;  (jai.i..,.  "  .UL,  C 

th>.:  as-.'  !  ;  ..  A  :  in  oquaticu.  ilC. 

S  -  .  ,  . .  :  .  .  ' 'u  i  i  ■  ■ !  I.-'.'  ;  i  A'l'ii  t . /  .  -  "c:,  ,j  i  .  .  '  n 

the  ■A  ■ ' '  ;  i  :  ;  roc  net  n  :s.  -  rd  1)0^1  a:,  e-  .,  1 1  i  .  n  5  a  n ; .  :  ■  ‘ 

'ill  i  .  :  .  ■  '  ; ■  • . .  i  g li'-'  1  t  or  ■■:!  i  ■’  i  u r'  I  n int'- 

t 1  ‘  i  '  A:  ,  n  r-  . 


,  n  r- 


■  .  .  ■  .  c.i-iluu  oy  L' "to  pC';'’iorn'i  a  special 

.^.iTcix  S'.u  - r  7  .v'nere  tlic-  leacinjj,  matrix  is  an  upper 
T;' i,.;i.,!. . .ij'  .atrix;  ciescriicU  ay  equation  '-b. 

;  :.a:i;'ou':  1  ii  ..  a’,  icc.  by  UbKEKi  to  pc,-i''li)rTr.  a  special  sym- 
'..aJCiC  ;.u'.  cipiy  .:a  neuooa  ix  equation  96. 

'i  -  ;  .ca.  .  0,0  .0.;.  oulli't;  by  EPLATE  to  t'>erform  a  lower  sym- 

i.ctri.,  .‘  o  o  ::,ui  t:. 'ly  oi  the  natural  .stiffness  components 
-ur:.".'  .:  .  Ti  Oaj.  integration,  as  in  equation  108, 

-  i  r". o  oa  -  lea  by  i'.'r.ATE  to  pei’:'orm  the  lower 

...itioi;  o:  tne  -tiifness  matx^ix  components 
O'n'i  '.r'  ;. . .1  in: egrat io.n ,  .-ts  describeci  by  equation  1C8 

-  a  s.t:  .  V.'.  ...a  oallca  by  CoLATE  to  determine  the  Gi.;irl- 
i.utici  c:  a;ea,  b;  noce,  for  the  quadrilateral  plate  elc- 
".e-.,.  i.nio  .:rea  function  it  needed  in  equations  112  and  1 

-  a  ;.  ,lr  .  on :  i.ie  called  by  CPLAi'E  to  iorm  a  transformation 
•.at:  ;:-;  fr  n:  nr.'.jr'a^  to  local  coordinates  as  described  in 


o.ibr  '  ;t  cai.Led  by  dPLATE  to  form  the  natural  to 
rlotal  c tor. i ' Hate  : ransf ormation  as  shown  in  equation  87. 
l:r.,::e  t.  x.'  foi^miat  ion  matrix  is  diagonal,  only  the 

'larorida  Sat  — m.-.ti'  >..  us  .;.re  x:torec. 

'  “  .i  .'  .1  r.'tU  :  .;.C  by  r'rL.blh.  to  '.'eX^l  om.'.  a  riatXta. 

',in..i.  : :  .m.  rmi.nAn  of  tl..'  .train-dLy;  i.xccrr.ent  a....  . 
ni  .q  n.-.i :  ;n  co.-.  ..scr'ibcd  in  equations  116  .inn  1  2 1  . 

-  .  C  i.;.i-eC.  l;v  ....-'.i'E ,  fOnl.T  ana  F  :''Moi !  1  i.i  - 

i  any  c*.  .’ci.e:  ..j  m.ati'ix  rnuj:  >-i'x  ..cation  o'  .nri'i  ^t  ra  nv  r.a:  .  c  , 
-•'i'Ctfd  :  I'or.  or  7>or..:.xon  to  ..iny  sub-mutrix  p... ait  ion. 

“  :  .J  n  raja  :  ..i.x:  C’ullral  aV  f  fh/.r.,  t  X'J  CUjC,;  I  utC  t  I'.O*  tait '  '. 

.  r'c,  .j.oa- ,  rciisr  inr  accelcr.-i t ion  anci  mass  vectc.rs  de.:.rj’'.  ;  . 
i  y  x:  ,aa  ;  i  a : a  .  .  ?  ,  119  ana  12  2. 

-  i  aabi’O'.tin.  c.;.  .ea  be  C 1  hATh  T'j  prepare  the  eii^iax;.;  .  : 
a.uj  .  .-atr  a;  q.ra-al  stit  inesa  matrix  with  elemont.,  :  rom  i:.i 
'r  .  :  iial  ■  1  ..ruJ  art  j  !  tor:.  L' .."lax  st:;  ■.nr:.:.  ::.a- 

T- :  i  .  in,.  ..>r  r>.  ■  •  s i ,k  ...ci'  i  bed  :  ,  alio:,  .  i.  1  . 


^  o  1  t he  ne  1  e ra  1  s  r  :  i  ;  ne  s  :na  t  r i.  x  i  ;i  t  o 
iaLai  L:rif  1  ne:.,.',  matrieen  .  'rhe  orip.inal 
mairix  ano  the  rina.!  t’xobal  ntifine;: 
:ne  came  mati’ix  space.  This  procedure  i 
:n,;tiL,n  103. 

i  ulled  by  i'TKANdf’  to  perrorm  a  i.pecia.i 
iroducr  used  !  i"!  stiifness  t  j 'ans  1  orma  t  ion 
:rfop:ris  an  e:  :  icienr  doul  matrix  m.:it  1 
^it>.;  of  t;.e  original  sub-matrix  as  aesc: 


called  by  PTRANSF  to  move  the  elements 
any  rank  and  ’,:'lace  them  into  new  matrix 


called  by  rTRAb'SF  to  move  elements  si:r. 
that,  the  receiving  sub-matrix  is  the  r.i 


dne  called  bv  CPLATE  to  perform  ,.n  cpecia 
such  that  the  post  multiplying  matrix 
;  shown  in  equation  126. 


■)K  :v 

RiriCATION 

irrAt-'  *  k  pr’o^i'cir.i  ^  nc;'.v  v^jr'  — 

:v;r  1  ■  lilv  includes  an  adnirLOiial  eli^ranr 
L  \ el  n^nt  )  '1  he  nri  i'.'i.nal  SAT  IV'  ; m:;. , 

..r  1,  iv.aoniJ  nnc  many  dilrei'ent  static  and 
: . .  n  i  V  s  2  a  .  ■  -  o  c  e  a  u  i"*  e  Cl ,  / 1  a  s  o  cn,  n  a  ±  t  e  a  d  v'  e  n  i  — 

.■  ;  he  new  elc-riont  ivas  chcclu  d  under  varxOuo 
conditions. 

v.ao  corr,par<“d  with  many  simpi-e  dege.nerate 
ems  ana  was  found  to  produce  excellent 
•ear.,  were  porformeu  for  element  TYPn  1  arid 
were  more  favorable  for  the  simple  cases. 

■.  .r..it.;  in  larger  displacements  while  TYPE  6 
.r\‘]er  than,  e.xact  solution  values, 
eet^on  aercribes  a  group  of  verif icat ion 
I .  r  / e  c  r  n e tit. d  a  d o u b  l y  c  rj r v e  c  blade.  . . .  -  e 
ontai.neu  bv  clansicc-.l  'plate  and  shell 
:a...tc  material  proper  t  i  eo  .  The  prO',  t.du:m 
:.ei'  be.' res  1C  aPi.>rox.m.-;Te  ..Oi  ;..t  1  o;._,  or  var lo 
cm  .a  1 1  lo.n  i- .  i.'nce  tni-  A,  .‘a  .r  n.-itri  ceSi  v.er 

s.'  n.I  j.  *  v,..i  j  S  d  j.‘  p  i  i  ^  V.  ^  ■.  1 .  .'.5  p  t.  - T*'v.^  A  -  i.4c;  w  i"!  .J.  .  u 

ri  .om’V’-t.  O':’  •'  L  r h I  - 1 ^-rrTt  vxp .  '/.csx 
c. :  r  . .  ■  I omo Ti t  .-cr  pt  ctiCn.  pc.i.c  v 

■■  at. .nice  .once  a  •.  -.c  tr.c  o..-..- 

: ecenae  i . 


4.  C 02/902 Jt  UNDER  UNIFORMLY  DISTRIBUTED  LATERAL  LOAD 
WITH  ALL  EDGES  SIMPLY-SUPPORTED 

(1)  Size  of  laminated  plate: 

a  =  10"  b  =  10"  t  =  4h  =  ,2" 


(2)  Properties  of  plate: 
p  =  .  0275^*’/in^ 


.  0275 

lb,.  3 
/in 

"3.3 

.  2 

.  2 

3.3 

0 

0 

■.135 

0 

0 

-.135 

L  0 

0 

r 

.67 

.  67 

11 

10^  lb 


1  Ox  10"  lb-in 
0  . 67j 


i'V 


(3)  Loading  condition: 
p(x,y)  =  p^ 

(4)  Boundary  conditions: 


X  =  0  ,  a ; 

w  =  0 

M 

II 

o 

V  = 

0 

X 

y  =  0,  b; 

II 

o 

M 

=  0 

u  = 

0 

y 

Deflection 

at  center. 

Wc; 

(W  ) 

=  114. 

4  X 

O 

1 

p 

in 

c  exact 

^o 

(w_) 

=  113. 

6  X 

10 

p_ 

in 

FIGURE  11 


[0/90]t  UNDi^R  FRLL;  VIBRATION  WITH  ALL  EDGES  SIMPLY-SUPPORTED 


(1>  Size  or  laminated  plate; 

a  =  10"  b  -  10"  t  =  2h  =  .2" 

(2)  Properties  of  plate: 
p  ^  .0275  Ib/in^ 


Ib/in 


( 5  ) 


1.65 

.1  Ol  6 

A  = 

.  1 

1,65  Ox  10 

.  0 

0  .  iJ 

3.37 

0  ol 

B  = 

0 

-3.37  0  X  10 

_  0 

0  oj 

13 . 8 

.  83  0] 

D  = 

.83 

13.8  Ox 

.  0 

0  . 8  3j 

(3) 

Loading  c 

ondition : 

free  vibration 

(4) 

Boundary 

conditions : 

lb 


X  =  0,  a; 
y  =  0,  b; 

P  requency : 
(f ) 


C  X  ci  C  L 


(f ) 


6.990  Hz . 
8.867  iiz . 


(5w  =  0 , 

6M  = 

X 

0, 

6v  = 

0, 

6N  = 

X 

0 

6w  =  0 , 

6M  = 

V 

0, 

6u  = 

0, 

6N  = 

V 

0 
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[90/0->/go]^  CURVKD  PLATE  UNDER  UNIFORM  PRESSURE  WITH  ALL 
L DG E S  ‘  S 1 M P LY - S U P PORTED 


(1)  Type  of  Blade:  (Figure  16) 

A  typical  blade  used  in  aircraft  engines,  named  J79  B/AL, 
was  approximated  using  64  quadrilateral  plate  elements  in  SAP4A 
(both  TYPE  6  and  TYPE  9).  An  equivalent  NASTRAN  model  was  also 
run  along  with  an  experimental  test  to  determine  fundamental 
frequency  at  zero  frequence  (results  from  AFSC  -  Wright-Patterson 
aFB)  . 

(2)  r'roperties  of  blade: 

.Material  used  corresponded  to  the  input  used  in  the  NASTRAN 
;^un  using  anisotropic  material  properties, 
heading  edge: 

s  =  .000407  lb  sec^/in  „ 

C:-;x  =  2  6.9E8  Ib/in^  Cxy  =  4.6E7  Ib/in'^ 

Tyy  =  20.3E8  Ib/in^  Gxy  =  7.3E7  Ib/in^ 

■Blade: 

p  =  .000251  lb  sec^/in 

oame  material  coefficients  as  above. 

(3)  Loading  Condition: 

Mass  and  Stiffness  distributions  for  eigensolution . 

(4)  Boundary  Conditions: 

Base  of  blade  completely  fixed  and  rest  of  blade  free. 

(j)  .Natural  Fr-equeiiC'/  of  Blade:  (Fiest  FLEX) 

(f)oxp  =  no  Hz. 

( !  )Na:'.tr  ..n  -  i  06  .  7  iiz  . 

(  f  )sap  (Tt.  )  -  1  08  .  3  Hz  . 

(f)sap  (T9)  =  L12.8  Hz. 


9 
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SECTION  V 

DISCUSSION  AND  CONCLUSIONS 


A  quadri  later.i j.  composite  plate  finite  element  has  been 
a -idea  to  the  SA?  IV  computer  prot^ram  library  to  be  used  on  plate 
tvpe  struotui^es.  The  element  is  a  four-noded  sub-parametric 
; iat  plate  element  excluding  transverse  shear  deformations.  The 
clement  is  "incompatible"  relative  to  mid-plane  surface  rota¬ 
tions  along  the  inter-element  boundaries.  The  element  converges 
rel.j:tively  well  despite  the  incompatibility,  as  long  as  the 
element  maintains  a  relatively  rectangular  shape  with  a  reasonable 
elerr.ent  aspect  ratio. 

The  element  was  run,  modelling  various  simple  plate/beam 
configurations  and  showed  excellent  results.  More  complex  models 
Ware  designed  to  test  the  composite  laminate  behavior  of  the 
.ioment,  as  described  in  section  IV.  The  models  included  flat, 
cylindrical  and  doubly-curved  shells.  The  results  obtained 
compared  favorably  with  classical  series  solutions,  (approximately 
1-6%  disagreement)  The  results  obtained  in  section  4.9  for  the 
typical  blade  conf iguration  shows  that  the  NASTRAN  and  SAP 
(TYPE  6)  element  values  were  slightly  lower  than  the  experimental 
number  while  the  result  of  the  new  element  (TYPE  9)  was  sligntly 
higher.  The  new  element  is  based  on  an  incompatible  formulation 
anc  in  general,  convergence  is  not  guaranteed.  But,  in  more 
c:ases,  the  elemenT  haS  been  found  to  be  slightly  stiffer  than 
compatible  elements  ^nd  therefore  the  frequency  is  higher. 

The  quadrilateral  clement  will  be  inherently  stiff  if  the 
foui'  points  do  not  represent  a  flat  surface.  Therefore  the  "quao" 
element  was  relaxed  to  better  represent  shell  behavior  by  allowing 
the  flat  plate  elemei.t  stiffness  coefficients  to  be  transformed 
relative  to  the  local  nodal  coordinates  of  the  element. 

The  thin  plate  theory  used  to  develop  the  element  does  not 
dc;coant  tor  normal  t'jrsional  effects.  Therelore,  non  global 
adiacent  elements,  when  assembled,  will  produce  a  singularity 


to  the  plate  it  the  two  elements  are  coplanar.  To  avoid 
ti;ia  singularixy  in  ^etieral,  an  artificial  torsional  stiffness 
or  scaffolding  matrix  was  added  to  each  of  the  four  nodes. 

This  does  violate  element  equilibrium  but,  if  the  magnitude  of 
tr.e  coefficients  are  maintained  relatively  "soft"  compared  to 
: he  plate  bending  characteristics,  overall  equilibrium  is  closely 
c..  I  -  n  t a  i  Tied  . 

Tr.e  composite  plate  element  can  be  used  in  the  existing 
-tatic  and  dynamic  analyses  contained  within  the  SAP  TV  program, 
it  ca,;.  be  used  witn  all  the  existing  elements  in  the  finite  ole- 
T.ont  library  as  long  as  the  model  effects  are  correct. 

The  element  was  not  being  developed  to  degenerate  to  a 
triangular  plate  element.  If  the  element  is  used  as  triangular, 
t;;c  local  effect  i:^  "too  stiff".  If  the  fourth  of  the  quacri- 
lateral  nodes  is  ixlaced  mid-plane  on  a  triangular  side,  a  better 
approximatic.n  can  be  obtained. 

The  compouirc  dement  presently  docs  not  contain  geometric 
tiiioning  cliccts.  ,.uch  as  those  required  in  high  speed  centri- 
:  ug.v.  nidcliinery. 

further  Devdoj .mont .. : 

F’urther  w.:>rk  :  i.ould  be  ..lone  Ct;  include  thu.  gt'ometric  s,t  i  : - 
•ar.ing,  c  1C  ;  ent '  io  account  ;  ot'  centripetal  accoleratioi- 
-d.'ecta  of  hign  .-p'dri  :)iade  systems.  Thir>  g.eometric  matrix  c/! 
coof  f  iciofits  would  allow  d  better  approximation  oi  blade  nencing 
stresses  and  closer  representat ion  of  blade  natural  frequencic: 
ut  high  spin. 

A  jire-j-rocessor  program  should  be  developed  to  hanc.li,  >;'i;mple 
uiaterial  laminates  as  a  function  o'  blade  position.  This  r;>  - 
gr'am  would  compute  the  A,  B  and  D  matrices  and  thermal  vector: 
needed  in  the  t;Af‘  program  and  produce  a  storage  file  for  ::tres;; 
recoverry.  This  program  should  also  plot  all  information  ica- 
iiisut  check  in, t’. 


A  should  be  Written  to  retrieve  istress  in- 

e.T.iri.'m  diid  .Tuite;' ial  information,  by  lamina,  to  be  used  with 
r- ;  .a  .T.. ; I  i on  output  and  curvatures  at  the  mid-plane  surfaces  to 
.u  iice  indivioual  l.imina  stress  to  be  used  in  failure  criteria. 
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appendices 


:  n.MX  A  -  IX-VT  TO  ELUILNT  TYPE  9  I:.  SA:  4A 

onij'/jot.o  [npi.it  to  tho  SAi-  IV  pronrani  dcncrin 
n.'W  element  9  defined  in  PA:  mA  ^ 

e  .jnu  .-iioulJ  De  appendotl  to  referenci;  1.  Tin- 
‘.'Yid,  9  i:.  nimilar  to  element  'I'YPE  6.  In  :  net 
j  replace  TYPE  f.  The  format  o;  t.  e 

n  is  eons.iGtent  with  the  SAP  IV  input.  Var;-.i 
nteper  must  be  right  -iustified  and  floating  po 
should  include  a  decimal. 


SV 


E  r-LATE  INPUT  TO  SAP4A 


TYPE  9  Composite  Pl.ite  Element  (QUADRILATERAL) 


:  .1  o  e  Colur.ir.3  Variable 

A.  Control  Cara  (615,  110) 


5  NPAR(l) 

6-lC  NPAR(2) 

11-15  NPARO) 

(1)  16-20  NPAR(4) 


C)  21- IS  NPAR(5) 


26-30  N?AR(6) 

(3)  31-40  NPAR(7) 


Remark 


Number  S 

Number  of  Plate  Elements 
Number  of  different  materials 
Material  Type  Key 
=0  Composite  Material  Prop. 

=1  Standard  Anisotropic  prop¬ 
erties  (sam.e  as  TYPE  6) 

Number  of  Global  Material  Vectors 
If  zero  or  blank  then  global  X 
direction  is  assumed  to  be 
material  x  axis. 

Integration  Order  (default  set 
to  2  ) 

Rotation  Stiffness  Factor 
(integer  number) 


b.  Material  Property  Information 


Two  types  of  material  can  be 
input  to  element  type  9;  ge.neral 
composite  material  and  aniso¬ 
tropic  material. 

Composite  Material  Properties  ( NPAR ( 4 ) . EQ . 0 ) 

Five  cards  must  be  input  for 
every  different  material.  (NPAR(3)) 


Cara  1: 

(110, 

20X, 

4F10. 

0) 

:  -10 

NN 

Material  identification  number 

11-30 

Blank 

31-40 

DEN 

Mass  density 

41-50 

AT(1) 

51-60 

AT(2) 

A„ 

thermal  vector 

components 

61-70 

AT(3) 

“  i 

(equation  66) 

Card  2  ; 

(6F10. 

0) 

1-10 

BT(1) 

11-20 

RT  ( 2  ) 

21-30 

BT  (  3  ) 

-  1 

thermal  vector 

component : , 

31-40 

DT(1) 

(equation  67) 

41-50 

DT(2) 

D„ 

thermal  vector 

components 

51-60 

DT(  3) 

•  i 

(equation  68) 

5fi 


C.ii'd  3  ; 


(6F10.0) 


A(l,l) 

il-20 

A(l,2) 

A 

21-30 

A(l,3) 

3 1  -  ‘4  0 

A(2,2) 

41-50 

A(2,3) 

51-60 

A(3 ,3) 

Card  4 : 

(6F10, 0) 

1-lC 

B(l,l) 

11-20 

B(l,2) 

B 

21-30 

B(l,3) 

31-40 

B(2 ,2) 

41-50 

B(2,3) 

51-60 

B(3,3) 

Card  5 : 

(6F10.0) 

1-10 

D(l,l) 

11-20 

D(l,2) 

D 

21-30 

D(l,3) 

31-40 

D(2,2) 

41-50 

D(2,3) 

51-60 

D(3,3) 

Matrix  coefficients  (upper 
triangular)  (equation  63) 


Matrix  coefficients  (upper 
triangular)  (equation  *64) 


Matrix  coefficients  (upper 
triangular)  (equation  65) 


3.2  Anisotropic  Material  Properties  (NPAR( 4 ) . EQ . 1 ) 

(4)  Two  cards  must  be  input  for  every 

different  material  (NPAR(3)) 


Card  1 : 

(110,  20X, 

4F10.0) 

1-10 

NN 

Material  identification  number 

11-30 

Blank 

31-40 

DEN 

Mass  density 

41-50 

AX 

Thermal  expansion  coefficient  ax 

51-60 

AY 

Thermal  expansion  coefficient  ay 

61-70 

AXY 

Thermal  expansion  coefficient  axy 

Card  2  : 

(6F10.0) 

1-10 

CXX 

Elasticity  element  Cxx 

11-20 

CXY 

Elasticity  element  Cxy 

21-30 

CXS 

Elasticity  element  Cxs 

31-40 

CYY 

Elasticity  element  Cyy 

41-50 

CYS 

Elasticity  element  Cys 

51-60 

GXY 

Elasticity  element  Gxy 

59 


Material 

Vectors  (25, 

5X,  3F10.0) 

NPAR(5) 

material  ve 

ctors  must  be 

input  (except  if 

zero  or  blank) 

1-5 

MV 

Material  vector  identification 

number 

6-lC 

Blank 

11-20 

DX 

X  direction 

cosine 

21-30 

DY 

Y  direction 

cosine 

31-^0 

DZ 

Z  direction 

cosine 

Element  Load  Multipliers  (5  cards) 
Card  1:  (4F10.0) 


1-10 

PA 

Distributed  lateral  load  multi¬ 
plier  for  load  case  A 

11-20 

PB 

Distributed  lateral  load  multi¬ 
plier  for  load  case  B 

21-30 

PC 

Distributed  lateral  load  multi- 
:;lier  for  load  case  C 

31-^0 

i'D 

Distributed  lateral  load  multi¬ 
plier  for  load  case  D 

Card  2  : 

(4F10.0) 

1-10 

TA 

Temperature  multiplier  for  load 
case  A 

11-20 

TB 

Temperature  multiplier  for  load 
case  B 

21-30 

TC 

Temperature  multiplier  for  load 
case  C 

31-40 

TD 

Temperature  multiplier  for  load 
case  D 

Card  3 : 

(4F10.0) 

1-10 

XA 

X-direction  acceleration  for 
load  case  A 

11-20 

XB 

X-direction  acceleration  for 
load  case  B 

21-30 

XC 

X-direction  acceleration  for 
load  case  C 

31-40 

XD 

X-direction  acceleration  for 
load  case  D 

60 


Care  ; 

( -1  Fi  0.0) 

1-10 

I'A 

Y-cii  re  cT  ion 
load  case  A 

acce- 1  i;rat  ion 

I  wi' 

1 

O 

Yn 

Y-n 1  recti  on 

accelerat ion 

r  e-:  ■ 

C.ii.C'  r- 

2  1  -  3  0 

YC 

Y-d i rect i on 

cicce  1  .o'd  t  ion 

^ '  y 

li-.ia  c.iSi-  i,‘ 

31-10 

e '  t 

V-cl  i  rect  1  or: 

I  Odd  case  D 

d  c  e  1  e  r'a  t  r  o  n 

:  cr 

Card  S  : 

(4ric. 0) 

1-10 

ZA 

7-direction 
load  case  A 

acce 1 era t iem 

fer 

11-20 

7,3 

7-direction 
load  case  E 

acce i erat ion 

for 

^o 

t-' 

1 

CO 

o 

ZC 

7-direction 
load  case  C 

acceleration 

:  o  r 

31-40 

ZD 

Z-direct ion 

acceleration 

f  or 

load  case  D 


: .  .  1 

leme^nt  Carets  ( 

315,  12, 

13,  12,  13,  15,  4F10.0) 

One  card 

for  each 

N  i  'AR  (2)  element. 

1-b 

NN 

Flement  number 

(5) 

f)-  1  0 

T 

Node  I 

11-13 

J 

Node  d 

T  L 

X  U  —  >d 

K 

Node  K 

7  ]  -  ?  S 

L 

Node  i. 

V  :  ) 

2  G  -  2  7 

NO; 

No.  of  Ociu, integration  :  o  1  it  t  . 

(7) 

2  S  -  3  3 

i'V 

Material  vector  ident  i  f  i  c  at  i  ,,n 
number 

(s) 

3  1  -  3  2 

Ii<;f.iSt: 

I’revious  Flement  r'e-use  coc<,' 

=  0  new  eleiTicnt 

=1  use  previous  element 

(9) 

3  3-3  S 

IM 

Material  identification,  number 

(13) 

36-40 

INCL 

Clement  generatictn  parameter 

(11  ) 

41-50 

TH 

Clement  thickness 

(12) 

51-60 

I’R 

i;iement  lateral  pressure 

(13) 

61-70 

TO 

Mean  temperature  variutior.  from 
the  reference  level  in  unoe- 
! ormed  pon i t i on . 

7  t-8  0 

TO 

Mean  1  entper'a  ture  nradicrit  ccros 

t  lie  ::)ie  11  I  ti  i  ekne:.:.  . 

N( .  t.e., : 

(1)  element  'i'YI'i,  ')  .illows  two  Mitlei'enl  I'e-rii  in  nit.  'iln: 

;  irst  form  is  tor  laminate  matrice;.  wlul'  tin  ..etaiml  i-.  I 'ii 
anisotropic  matric'cs.  The  Later  tor'm  i:  idetirical  to  cle¬ 
ment  TYn;  6  int>ut. 


t.l 


(.'2)  A  rual.iM'i.il  axi::  mur.t  1m?  dc-'fined  relative  to  the 

material  propar  t  i  ea  f  ornuHl  in  t  lie  A,  B  and  D  nuitric'es. 

(3)  KiMtational  Stiltnesn  Factor  i r.et  by  multiplyinp,  NI'ARCV) 
times  l.F-8.  befault  foi'  NlAk(7)  if,  100. 

(4)  Material  input  in  this  election  is  identical  to  thcjt  of 
element  I’YPL  0. 

(5)  The  T,J,K  and  I,  indices  dcline  the  element  connectivity 

iiul  also  the  element  normal.  The  element  "z"  coordinate 
is  !  or'med  by  the  rip.ht  hand  rult>  as  T  p.oer,  to  J  p.oes  to  K, 
etL'.  The  element  local  ,ixi:'.  is  determine^l  by  the  [-.r'oiec- 
ti<.)n  oi  the  I’lobal  materi.il  axis  onto  the  element.  Once 
the  x-ax  is  is  'ietermined,  the  local  y  if,  fot'med  from  z 
■ind  X.  All  stress  output  is  in  this  relerc'nce.  If  node  L 
e'jual  K  or  il  zero  or  let  t  bl.ink,  the  pt'opi'am  will  assume 
tliat  the  clement  is  tri  aiu’.u  1  ar .  The  resultinp  local  stiff¬ 
ness  is  then  "too  stiff". 

( i'  )  The  number  ot  <',auss  intep.ra  t  i  on  points  can  vary  .is  2  or  3. 

It  a  v.ilue  is  -.et  abctve  or  below  these  numbers,  the  f?rop,r'dm 

will  reset  it  to  ?.  The  delault  valut'  is  set  FFARCG). 

(7)  The  f.lubal  material  vcc'tor  must  be  p.roatLU'  th.tn  oi'  eipi.tl 

to  1  aii'l  less,  than  or  e<)u.il  to  N1'AK(5).  II  NFAKfl))  is  bl.ink 
OI'  zert.-),  then  NI'AK(li)  is  set  e<|ual  to  1  .ind  the  plobal  vec¬ 
tor  is,  .ilipned  .ilont'  the  I’h'b.il  X  .ixis..  befault  value  it; 
set  to  1 . 

(3)  if  rin  el.meiit  h.is.  the  saiiK'  i'i.in.ir  t.ize,  s,,iine  ot' i  en  1  .t  t  i  c'li  in 

.S[',ice  ani  the  s.ime  oli'mt'nl  l.'.idinp,  ixii'.imc't  ers  as.  the  pre¬ 

vious.  ele;;i#.'ni  ,  then  settini’  I  KFUF.F  oqu.il  to  ]  will  usi;  tlu' 
s,.ime  plob.il  eli'mcnt  ;;t  i  f  i  ne:'.s.  ,ind  load  vector's  for  ass(,'nib]v. 
i'efault  is  .set  to  n. 

(3)  The  mater  i-il  lb  numbitr  muf.t  l.e  betw.'cn  1  .ind  NI'AFfl).  De- 
lault  is  set  to  1. 

(in)  idemeiit  flenerri  t  i  on  p.tr,  ime  ter ;  Flemeni  c.irds  must  be  in 

'■lernent  number'  s,(’quence.  It  .'lement  r'.irils  .jri'  omitted,  the 
pt'opr’am  wi  I  1  cenerate  the  mis.siii)',  cai'ds,  ,is  lot  lows: 

The  increment  lor'  the  element  number  i:  one. 

i  .  ,  bh  ■  ,  ,  -  hli  •  ■*  1 

1  M  I 

Tip,  t  ■  r ,  r '  t 'e  ■ ,  I  ■!  rir  I  1  I  n’  iin'romi-nt  I'.r  ni>d.il  eonne.l  i  v  i  I  is  IN't'I,. 


‘i  +  1 


+  INCL 


J.  + 

1 

INCL 

K.  + 

1 

INCL 

4.1  ^ 

L.  + 

1 

INCL 

If  INCL  is  left  blank 

then 

INCL  is 

set  to  1.  Material 

identification,  elemen 

t  thi 

ckness , 

distriDuted  lateral 

load,  temperature  and  temperature  gradient  for  the  element 
are  then  the  sam.e  as  for  the  first  element  in  the  generatea 
group.  The  last  element  card  must  be  input  to  exit  element 
group  properly. 

(11)  The  plate  thickness  is  used  mainly  to  compute  the  mass  of 
the  element.  Default  is  set  to  1.0. 

(12)  The  pressure  is  normal  to  the  surface  of  the  element.  The 
positive  direction  for  the  pressure  loading  vector  is  in 
the  positive  direction  of  the  local  z  coordinate. 

(13)  The  tem.perature  required  is  the  mean  temperature  difference 
(TO)  from  the  reference  temperature  of  the  element  in  a 
undeformed  state.  The  TG  is  the  mean  thermal  gradient 
through  the  element  thickness. 
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APi-ENDIX  B 


INPUT  TO  rRL’-PROCL'SSnR  F'ROC.kAK  LAYUP 


A  i^re-proc.  :>:'.L)r  prop,r>-im  wa:,  Uc‘Veloi)ed  to  calculate  the 
,  P.  and  D  rnati’ici-..  j.'.  tieucribeU  in  equalionu  63  thcouph  63  ai.d 
t  e  rh<.r'r:i>il  loac  v.n.  !  oru  .Icucribcc  in  oquaticn;.  66  tJirfiup.ji  68. 

:ie  iT.aterial  Ln ;  crr:iat  i  on  by  layei'  i  r.  input  r.c  iamina  -position 
-i.'id.  :  iber  rn^ientat  ion .  Tne  pre-processor  i  hca  computes  the;  C 
r: rices  involvl:.,  the  tensor  transformation  and  coiistructs  tl'ie 
eioment  materia’,  m^utrix  Li  rel.jtive  to  the  mid-plane  of  the 
plate.  Similar  the;rmal  vectors  are-;  also  computed  and  oufput  . 
Tnis  information  ’is  then  used  directly  for  SAPhA  -  TYPE  9. 


flUAK-woT 


IK'FUT  TO  PROGRAM  LAYUP 


A. 

Number  of  Cases 

Card  (12) 

: .  o  1. 1 

Columns 

Variable 

Remark 

{!) 

1-2 

N CASES 

Enter  the  total  number  of 
laminate  configurations  to  be 
considered . 

B . 

Heading  or  Titl 

e  Card  (8A10) 

(2) 

1-80 

ITITLE 

Enter  the  title  information  to 
be  printed  with  the  output. 

C. 

Laminate  Data  Card  (415) 

(3) 

1-0 

NLAM 

Enter  the  number  of  laminae  in 
the  laminate. 

(‘t) 

o-lC 

NMAT 

Enter  the  number  of  different 
materials  in  the  laminate. 

11-15 

I  TRAN 

If  ITRAN  is  zero  or  blank  the 
transformed  thermal  properties 
of  each  lamina  will  not  be  part 
of  the  output. 

16-20 

lORD 

If  lORD  is  left  blank  the  lami¬ 
nate  ordinates  will  not  be  output 

J . 

y.aterial  Proper 

ty  Card  (6F10. 

3) 

(b) 

1-10 

ElKK) 

Enter  lamina  modulus  in  fiber 
direction 

••  1  _  ) 

i;22(K) 

Enter  iamini  modulus,  in  direction 
transverse  Lo  fibers. 

: .  -  .5  0 

XNKK) 

Enter  manor  Poisson's  ratio. 

-  4  0 

fill  (K) 

Enter  lamina  ;;.hear  moculus. 

-  50 

THERMKK) 

Enter  C.T.;,  .  in  ;  iber  direction. 

51-60 

THERM2(K) 

Enter  C.T.E.  in  ;  ransver.oe 
direction . 

L . 

uamina  Data  e-ir 

d  (no. 3,  no. 

PIO. 3) 

1-  .  0 

T(J) 

Enter  thickness  of  the  lamina 

11-20 

MATL(  J) 

Enter  the  nvimher  that  ^^dent  i f  ics 
the  material  of  the  ‘  lamina^^_ 

Ent(;r  the  orientation  of  the  J 
lamina  with  resj>ect  to  the 
laminate  axe:;. 

2  1-30 

FJKJ) 

Notes : 

(1)  There  are  no  program  restrictions  on  the  number  of  cases 
that  may  be  analyzed  in  a  single  run. 

(2)  Begin  each  new  data  case  with  a  heading  card. 

(3)  The  program  is  currently  capable  of  handling  up  to  48  lami¬ 
nae  per  layup.  This  can  be  increased  by  changing  the 
appropriate  dimension  statements  as  shown  in  the  program 
listing. 

(4)  Although  no  limitation  on  the  number  of  different  materials 
need  be  imposed,  the  program  dimension  statements  currently 
allow  for  a  maximum  of  NMAT  =  4.  However,  this  can  also  be 
increased  if  necessary. 

(5)  This  material  data  is  vendor  information.  One  card  is 
required  for  each  material  (i.e.,  K  =  1,  NMAT). 

(6)  One  card  is  required,  per  layer,  in  the  laminate.  (i.e., 

J  =  1,  NLAlf). 


APPENDIX  C 


COMPUTER  RUN 


INPUT  AND  OUTPUT 


This  appendix  contciins  the  complete  input  and  resultinp, 
output  of  the  SAP'^A  proj^ram  using  the  new  composite  plate  e  K  - 
ment  (TYPE  9).  The  example  is  the  model  contained  in  section 
4.7  (Curved  Plate  under  Uniform  Pressure).  The  first  secti.  ii 
contains  a  listing  of  the  card  images  used  to  execute  the 
gram.  The  second  section  is  the  complete  output  resulting 
from  this  input. 
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EQUATION  NUMBERS 
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